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Abstract 



We have observed the cold collision frequency shift of the 1S-2S transition 
in trapped spin-polarized atomic hydrogen. We find Ai>is-2S = —3.8 ± 0.8 x 
10~ 10 n Hz cm , where n is the sample density. From this we derive the 15- 
2S s-wave triplet scattering length, a\s-2S = —1-4 =L 0.3 nm, which is in fair 
agreement with a recent calculation. The shift provides a valuable probe of 
the distribution of densities in a trapped sample. 
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Methods for cooling and trapping neutral atoms have opened a low energy regime where 
collisional phenomena involve only a single partial wave. The ground state s-wave scattering 
lengths for many of the alkali metal atoms have been obtained from combinations of theory 
and experiment. Those studies were motivated in part by the need to estimate evaporative 
cooling rates H, study the stability of Bose-Einstein condensates 0|§, and understand 
collisional shifts of atomic fountain frequency standards [|].|| . 

We have measured a density-dependent frequency shift of the two-photon 1S-2S tran- 
sition in trapped atomic hydrogen. This allows us to deduce the 1S-2S triplet scattering 
length, which has been calculated by Jamieson, Dalgarno, and Doyle 0. The frequency 
shift also provides a probe of the density distribution of our trapped sample, and constitutes 
a valuable tool for the study of Bose-Einstein condensation. (See the accompanying paper 
on BEC in atomic hydrogen j7|.) 

Interactions between neighboring atoms shift and broaden atomic spectral lines. For 
bosons at low temperatures, where only s-wave collisions occur (below about 1 K for H ||), 
the elastic contribution to the shift can be described through the density-dependent mean 
field level shift @ 

8nh 2 an , . 

E = . 1 

m 

Here, a is the atom's s-wave scattering length for collisions in the gas, and all particles 
have mass m. The gas has density n, and is assumed to be nondegenerate. The difference 
between the mean field energy in the excited and ground states results in a frequency shift. 
For weak excitation of the 1S-2S transition in spin-polarized hydrogen, 

Az/i5„2S = («15-2S — dlS-ls) — , (2) 

m 

where the scattering lengths are for the triplet potentials. 

Additional contributions to the shift can result from inelastic processes. The most im- 
portant corrections should arise from collisional hyperfine transitions and quenching of the 
25* state. Hyperfine transitions in the spin-polarized sample arise only through weak mag- 
netic dipole interactions JTO], and in our experiment we observe no evidence for collisional 



quenching on a millisecond time scale. Thus it is reasonable to neglect inelastic processes 
when describing the observed density-dependent frequency shift. 

Analogous frequency shifts are observed in hydrogen masers JTl| and cesium atomic 



fountains |[L2| . Using the collision-based formalism commonly applied to these systems, one 
can derive Eq. 2 and also show that homogeneous broadening due to elastic collisions is 
negligible in the present study because the phase shift per collision, 27ra/Ad e Brogiie5 is small. 



Our techniques for cooling and trapping hydrogen have been described elsewhere (13 
and are only briefly summarized here. Molecules are dissociated in a cryogenic discharge, 
and the atoms thermalize by collisions with the 250 mK liquid 4 He-coated cell wall and each 
other. A fraction of the "low-field seeking" F = l,mp = 1 atoms settle into the magnetic 
field minimum of a 500 mK deep Ioffe-Pritchard magnetic trap \TM (Fig. 1). An axial bias 



field of a few gauss in the trapping region inhibits nonadiabatic spin flips. 

The wall temperature is decreased to 120 mK. This increases the wall residence time of H 
atoms, facilitating recombination, so that contact with the wall effectively removes energetic 
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atoms. This process initiates evaporative cooling (T^JT^| and thermally disconnects the 



sample from the wall. The evaporation is forced by decreasing the trap depth, which is set 
by a saddlepoint in the magnetic field. Sample temperatures as low as 110 /iK and densities 
approaching 10 14 cm -3 are achieved. Under these conditions the sample is a needle-shaped 
cloud with a diameter of 200 /im and a length of 8 cm. 

The speed with which the trap threshold can be lowered is limited by the need to main- 
tain thermal equilibrium through elastic collisions during the evaporation. After the forced 
evaporation is terminated, the sample temperature remains constant. Its value is deter- 
mined by the balance between cooling due to evaporation of atoms over the threshold and 
heating due to dipolar spin-flip collisions. Dipolar decay preferentially removes atoms from 
the bottom of the trap where the density is highest. 

The temperature |T7| and density [16j are measured by rapidly lowering the trap threshold 



to zero while monitoring the power deposited on a bolometer. The power results from 
molecular recombination and is proportional to the flux of atoms escaping from the trap. The 
dependence of the flux on the threshold reveals the energy distribution of atoms in the trap, 
which, combined with a knowledge of the trapping fields, yields the sample temperature. 

We measure the density by releasing the atoms from the trap after holding them for 
different times t following the forced evaporation. The total recombination energy deposited 
on the bolometer is proportional to N(t), the number of atoms remaining in the trap at time 
t. N(t) decreases with time due to dipolar decay, which removes atoms from the sample at a 
density-dependent rate (Fig. 2). This process obeys the two-body equation dn/dt = —gn 2 , 
where n is the local density and g = 1.1 X 10~ 15 cm 3 /s is the calculated rate constant |T0| , p!8| . 
When averaged over the density distribution in the trap, this implies that N(t) obeys 



N(0)/N(t) = 1 + ngn (0)t. (3) 

Here, no(t) is the maximum density in the sample at time t, and k is a numerical factor which 
depends on the trap geometry. Computations show that k is typically in the range 0.2-0.3. 
For a given geometry, the uncertainty in k is 10%, arising from imperfect knowledge of the 
magnetic fields. The number of atoms for a given trapping configuration is reproducible to 
a few percent. 

Once the sample temperature and density are known, we measure the cold collision 
frequency shift by two-photon Doppler-free 1S-2S spectroscopy [y|. The 243 nm laser that 
excites atoms to the metastable 25 state is stabilized by an optical cavity and has a linewidth 
of about 1 kHz. 8 mW of radiation is focused to a waist radius of Wq = 50 fim in the trap. 
The beam lies along the axis of the atom cloud and is retro-reflected by a spherical mirror 
to provide the standing wave required for Doppler-free excitation. A mechanical chopper 
pulses the laser beam at 1 kHz with a 50% duty cycle. Following each pulse, an electric field 
of 10 V/cm Stark-quenches the excited atoms by mixing the 2S and 2P wavefunctions. This 
causes prompt radiative decay by emission of a Lyman-a photon (122 nm). The emitted 
photons are counted using a microchannel plate detector located at the end of the trap, 
behind the retro-reflecting mirror. Due to the low collection solid angle (2 x 10~ 2 ster), 
absorptive losses, and detector quantum efficiency, the total detection efficiency is only 
10~ 5 . 

Typical data for determining the cold collision frequency shift is shown in Fig. 3a. The 
density is highest immediately after the forced evaporation, so the first recorded spectrum 
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is shifted the most. Due to the distribution of densities in the sample, the line is in- 
homogeneously broadened. The density decreases on a 20 s time scale primarily due to 
collisions with helium atoms evaporated from the mirror surface by the laser; hence suc- 
cessive scans probe different densities. The later spectra exhibit the exponential lineshape, 
exp(— | z/ 1 /5z/ transit) f2(J, expected for Doppler-free two-photon excitation by a Gaussian laser 
beam of a low density sample. Here, v is the laser detuning from resonance at 243 nm in 
Hz, bvt r ans%t = y ^ksT / 'm j ' (4ttwq) is the linewidth due to the finite interaction time of an 
atom with the laser beam - typically 2-5 kHz. Measurement of the low density linewidth 
provides an independent check of our sample temperature . 



To quantify the cold collision frequency shift, one must account for the distribution of 
densities in the trap. Numerical simulations of the spectra show that the shift of the line 
center is about 0.75 of the shift associated with the maximum density in the sample, n^. 
This factor depends on the trapping magnetic fields, sample temperature, and laser focus 
position with respect to the atom cloud. There is a 10% uncertainty in this correction for a 
given configuration. 

The integral of the 1S-2S spectrum is proportional to no(t), and it decreases smoothly 
in time. This allows us to extrapolate from the initial density before the first scan, found 
as described above, to determine the density for each successive scan. A fit to data such as 
shown in Fig. 3b yields x-, where x = &- v is-2s/ n - Measurement errors are small compared 
to the mentioned systematic uncertainties. The long term drift of the laser frequency, ~ 10 
kHz/hour, is taken into account. Note that no knowledge of absolute frequency is required 
to determine \ from a given series of scans. 

To examine the sensitivity to systematic effects, x was measured in different trap config- 
urations, with sample temperatures between 110 and 500 fiK and initial maximum densities 
in the range (2 ~ 7) x 10 13 cm -3 . In Fig. 4, the various measurements are plotted versus 
sample temperature. We observe no significant temperature dependence of the shift, consis- 
tent with theory. Any significant nonlinear density dependence of the shift would manifest 
itself in data such as Fig. 3b, but none is evident. 

From a weighted average of the various measurements, we find x = —3.8 ± 0.8 x 
10~ 10 Hz cm 3 . In an apparatus optimized for spectroscopy with better known magnetic 
field and laser geometries, the uncertainty could be greatly reduced. A revision of the cal- 
culated value for g would cause a proportional change in the value of x- Experiments have 
verified this decay constant at the 20% level pi |, but since the calculation involves well 



known ground state potentials, the calculated value is expected to be more certain. 

Equation 2 relates x to the s-wave triplet scattering lengths, assuming inelastic processes 
can be neglected. From this, we derive ais~2S — — 1.4±0.3 nm. We have used the theoretical 



value of ais-is = 0.0648 nm p2|, which constitutes only a small contribution to x- The 
value calculated by Jamieson et al. ||, a\s-2S — —2.3 nm, is in fair agreement with this 
measurement. 

Having calibrated the cold collision frequency shift, we can use two-photon spectroscopy 
to measure sample density provided the shift is comparable to or greater than the transit- 
time width, (^transit- The spectrum, 5(z/), is given by the approximate expression 

3(u) oc J d 3 r n(r)I 2 (r)5(u - xn(r)), (4) 

which neglects atomic motion and laser linewidth. J(r) is the laser intensity, and the delta 



4 



function enforces the resonance condition. In our experiment, at high densities and low 
temperatures, and especially in the regime of Bose-Einstein condensation, this probe is more 
sensitive and accurate than the bolometer. In addition, it provides valuable information on 
the density distribution |]]. 

1S-2S spectroscopy of cold trapped hydrogen can be used for optical frequency metrology 
and precision measurements. Linewidths of about 1 kHz have been achieved in an atomic 
beam |]23|| , but with the long coherence time possible in a trap [|B| it should be feasible to 
approach the 1.3 Hz natural linewidth of the transition. To reduce the effects of the cold 
collision shift below this level, it is desirable to work at densities of < 10 10 cm -3 . Even with 
such a decreased density, large signal rates are still possible in an apparatus optimized for 
optical access. 
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FIGURES 




L 243 nm laser beam 

FIG. 1. Schematic diagram of the apparatus. Coils create a magnetic field with a minimum 
along the trap axis, which confines the sample. The 243 nm laser beam is focused to a 50 /^m 
beam radius and retroreflected. A 500 /j,s laser pulse promotes some atoms to the metastable 2S 
state. An electric field then Stark-quenches the 25 atoms, and the resulting Lyman-a fluorescence 
photons are counted by the microchannel plate. Not shown is the trapping cell which surrounds 
the sample and is thermally anchored to a dilution refrigerator. The actual trap is longer and 
narrower than indicated in the diagram. 
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FIG. 2. Dipolar decay of trapped hydrogen. The sample density is found from the slope of 
N(0)/N(t), the inverse of the normalized total number of atoms remaining in the trap. The data 
shown indicates a density of 6.0 x 10 13 cm~ 3 . 
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FIG. 3. (a) Cold collision frequency shift observed in the spectra of a single 120 /iK sample 
with initial maximum density of n-o(O) = 6.6 x 10 13 cm~ 3 . The furthest red-shifted spectrum was 
recorded first, and subsequent spectra were taken at decreasing sample densities. Each data point 
represents 12 ms of laser excitation and one scan is recorded in 0.5 seconds. 20 seconds elapsed 
between the first and last scan shown. The smooth curves are the result of numerical simulation. 
Only 5 of the 40 spectra taken of this sample are shown. 

(b) Frequency shifts from spectra of sample described in (a) after correcting for inhomogeneous 
sample density. For this data, the linear fit yields x = —3.30 ± 0.6 x 10~ 10 Hz cm 3 , where 
X = Az/i5_2s/n. (Note, Avisos, the shift in the transition frequency, is twice the shift in the 
frequency of the 243 nm laser.) 
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FIG. 4. The frequency shift parameter x, determined as described in Fig. 3 for various trap 
configurations. No significant dependence on sample temperature is observed. The error bars 
reflect systematic uncertainties in the magnetic trapping fields and laser geometry, with small 
contributions from measurement error. The dashed line is the weighted mean of all measurements, 
X = —3.8 ± 0.8 x 10" 10 Hz cm 3 , and the quoted uncertainty is indicated by the double arrows. 
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